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Abstract 

A new phenomenological cluster-hadronization model is presented. Its specific fea- 
tures are the incorporation of soft colour reconnection, a more general treatment 
of diquarks including their spin and giving rise to clusters with baryonic quantum 
numbers, and a dynamic separation of the regimes of clusters and hadrons according 
to their masses and flavours. The distinction between the two regions automatically 
leads to different cluster decay and transformation modes. Additionally, these as- 
pects require an extension of individual cluster-decay channels that were available 
in previous versions of such models. 
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1 Introductory note 



Multi-hadron and jet production in high-energy particle reactions is a ba- 
sic property of the strong interaction [1,2]. A successful description relies on 
a factorization, which permits the separation of the perturbative evolution 
from the non-perturbative development of an event. The perturbative regime 
can be characterized through calculations of hard matrix elements and sub- 
sequent multiple parton emissions - the physically appealing parton-shower 
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picture 1 . However the entire hadron-production mechanism cannot be pre- 
cisely predicted because of the lack of the understanding of non-perturbative 
QCD effects, i.e. hadronization. For the transition of a coloured partonic sys- 
tem into colourless primary hadrons, this implies a need for phenomenological 
models. Lastly, after the primary-hadron genesis, decays of unstable hadrons 
are accomplished. Employing the separation ansatz, Monte Carlo event gener- 
ators such as JETSET/PYTHIA [5] or HERWIG [6] proved to be a successful tool 
for the description of multiparticle generation in high-energy physics. 
Concerning the transition process, such Monte Carlo schemes are either based 
on the Feynman-Field or independent fragmentation [7], on the Lund string 
[8] and UCLA [9] model (JETSET/PYTHIA), or on the cluster-hadronization 
model (HERWIG). The latter concept 2 , initially proposed by Wolfram and Field 
[11,12], and further advanced, among others [13], by Webber and Marchesini 
[14,15], explicitly rests upon the preconfmement property of QCD [16] and the 
LPHD hypothesis [17]. Such cluster models are usually formulated in terms 
of two phases: cluster formation accomplished through the non-perturbative 
splitting of gluons left by the parton shower into quark-antiquark pairs, and 
cluster decays leading to the additional creation of light-flavour pairs. 
To understand the physics at present and future colliders, e.g. the Tevatron 
at Fermilab and the LHC at CERN, one fundamental cornerstone is the im- 
plementation of new Monte Carlo event generators, e.g. PYTHIA7 [18-20], and 
HERWIG++ [18,21,22]. The development of the C++ Monte Carlo event generator 
SHERPA (Simulation of High Energy Reactions of PArticles) [18,23] is a step 
in the same direction. The modified phenomenological cluster-hadronization 
model presented in this paper contributes as a further module to the construc- 
tion of the SHERPA package. The basic features of the new model are: 
Soft colour reconnection is accounted for in the formation and decay of clus- 
ters. The flavour-dependent separation of the cluster regime from the region 
of hadron resonances yields the selection of specific cluster-transition modes. 
The two regimes are distinguished by comparing the mass of the cluster with 
the masses of the accessible hadrons matching the cluster's flavour structure. 
So far, the cluster scheme presented here is implemented only for electron- 
positron annihilation, and, for simplicity, only the light-quark sector is con- 
sidered. An extension to heavy quarks, however, is straightforward. 
The paper describing our cluster-hadronization model is organized as follows: 
first, different aspects of cluster formation are discussed in Sec. 2. Subse- 
quently, in Sec. 3, the parametrization of light-flavour pair creation is pre- 
sented. The model's description is concluded by exhibiting cluster transfor- 
mation and fragmentation processes, which lead to the emergence of primary 
hadrons, see Sec. 4. The first results obtained with the new hadronization 
scheme are shown in Sec. 5 for the process e + e~ — > 7*/Z° — > dd,uu,ss — > 



Perturbative QCD cascades can be formulated in two complementary ways, either 
in terms of quarks and gluons or in terms of colour dipoles [3,4]. 
2 Recent developments may be found, e.g. in [10]. 
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hadron jets. 



2 Cluster formation 

The parton shower describes multiple parton emission in a probabilistic fash- 
ion [2]. By factorizing the full radiation pattern into individual emissions it 
employs the large- Nq limit of QCD. This organizes a binary tree, i.e. a planar 
structure, of the partons. It also ensures that, once the colour structure of the 
initial partons from the hard matrix element is fixed, the colour structure of 
the partons at the end of the parton shower is unambiguously determined. 
In our model, the non-perturbative transition of these partons into primary 
hadronic matter, clusters, is accomplished by the following steps: 

(1) To guarantee the independence of the hadronization model from the 
quark masses eventually used in the parton shower and to account for 
a gluon mass needed by the model, all partons are brought to their con- 
stituent masses [14], 0(0.3 GeV), 0(0.3 GeV) and 0(0.45 GeV) for u, d 
and s flavours, and 0(1 GeV) for the gluon, respectively. For this tran- 
sition a numerical method, involving several particles and consisting of a 
series of boosts and scaling transformations, is employed. However, these 
manipulations are applied only to parton-shower subsets that are in a 
colour-singlet state. 

(2) Since in cluster-hadronization models the clusters consist of two con- 
stituents in a colour-neutral state made up of a triplet-antitriplet, the 
gluons from the parton shower must split (at least) into quark-ant iquark 
pairs [12]. So, a transition - in principle non-perturbative - transition 
g — > qq, DD into a light quark-antiquark pair qq or a light antidiquark- 
diquark pair DD (see Sec. 3) is enforced for each gluon. The respective 
flavour composition of the gluon's decay products is obtained with the 
same mechanism as used for cluster decays; see Sec. 3. Quarks or di- 
quarks that cannot be produced owing to too high masses are discarded. 
The kinematical distribution obeys axial symmetry; the energy fraction 
z of the quark (antidiquark) w.r.t. the gluon is given by a density propor- 
tional to z 2 + (1 — z) 2 , i.e. the gluon splitting function 3 . The limits on z 
are fixed only after the flavour of the decay products has been selected. 

(3) In contrast to the Webber model of cluster fragmentation [24], our model 
may also incorporate soft colour reconnection 4 effects by eventually re- 
arranging the colours of the partons forming the clusters. Starting with a 
simple cascade, Fig. 1 schematically shows the two options to arrange two 

3 Obviously, for antidiquark-diquark pairs, this is a simplistic assumption, since it 
neglects, at least, the different spin structure of diquark production. 

4 Other soft colour reconnection models are presented, e.g. in [25,26]. 
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Fig. 1. Both options of cluster formation for a minimal qgq -^qq'q'q cascade. The 
zig-zag lines connecting the quark lines symbolize the soft exchange of colour quan- 
tum numbers, which is responsible for the colour reconnection. 

colour neutral clusters out of four quarks or diquarks. The first - direct 
- case corresponds to the usual cluster formation and reflects the leading 
term in the 1/iVc expansion. The second - crossed - configuration keeps 
track of subleading terms. Motivated by the well-known colour suppres- 
sion of non-planar diagrams w.r.t. planar ones, the relative suppression 
factor due to colours is taken to be 1/Nq. Additionally, a kinematical 
weight is applied for each of the two possible cluster pairings. For the 
pairing ij, kl this weight reads 



where the quantity t , of the order of 1 GeV 2 , denotes the scale where the 
parton-shower evolution stops and hadronization sets in. As a measure, 
Wn functions such as the invariant mass 



of the parton pair (and therefore the cluster), or their relative transverse 
momentum, similar to the Durham k± jet-scheme [27]: 



might be used. The actual colour configuration of the considered system 
is then chosen according to the combined colour and kinematical weight. 
Ultimately, this reshuffling is iteratively applied to combinations of two 
colour-singlet pairs of partons in the colour-ordered chain. 
Of course, users who are not interested in colour reconnection have the 
possibility to switch this option off. 
(4) The cluster formation is accomplished by merging two colour-connected 
partons, quark or antidiquark and antiquark or diquark, into a colourless 
cluster. In this way, four different types may arise, mesonic (qilfr and 
DD'), baryonic (qiD), and antibaryonic (Dq^) clusters. The total four- 
momentum of these clusters is just given by the sum of their constituent 
four-momenta [12]. 







(2) 




(3) 
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3 Parametrization of light-flavour pair production 

In our model the gluon splitting at the beginning of the cluster formation 
phase and all cluster decays rely on the emergence of light-flavour pairs, see 
[12,14]. During hadronization, which typically sets in at a scale of 1.0 GeV, 
there is no possibility for heavy-flavour pair generation [8]. The appearance 
of baryonic structures is tied to the creation of light diquark-antidiquark 
pairs 5 . In contrast to the Webber model, in our approach the total diquark 
spin S is explicitly considered. Thus, qq and D$Ds, where q E {d, u, s} and 
D s G {ddi, ud , udi, uui, sd , sdi, su , sui, ssi}, occur as the possible pairs. 
Apart from their masses influencing their emergence, the created pair func- 
tions only as a flavour label. Furthermore, the pair generation is assumed to 
factorize, i.e. to be independent of the initial flavour configuration. There- 
fore, the only interest lies in finding suitable pair-production probabilities, 
i.e. flavour and spin symmetries should be correctly respected and reasonable 
hadron multiplicities should be finally obtained in the hadron production. 
In our model a phenomenological parametrization is achieved by employing 
hypotheses leading to a general "flavour dicing" scheme. This scheme is ap- 
plied to both regimes, cluster formation and decay. The hypotheses are: 

(1) The emergence of diquarks, i.e. baryons, is suppressed through a factor 
p B with < ps < 1. 

(2) SU(3)^ symmetry is applied, but is assumed to be broken. This is mod- 
elled by a strangeness suppression parameter p s , < p s < 1, whereas 
the production of u and d flavours is equally probable (strong-isospin 
symmetry), hence p u ^ = (1 — p s )/2. As mentioned above, p c ^ = 0. 

(3) Spin and flavour weights: the spin-S 1 diquark states (S = 0, 1) get a weight 
proportional to 25 + 1. Additionally, a combinatorial factor of 2 and 1 is 
applied, depending on whether different or equal flavours constitute the 
diquark. But the fact that all states in the baryonic SU(3)^ octet and 
decuplet appear equally likely has to be reproduced. This gives rise to 
extra weights on the individual diquark types. In particular, the combined 
diquark weights read (up to the baryon suppression factor): 

w D=ud,sd,su = VDi w D=ud,sd,su = %PD, w D=dd,uu,ss = ^PDi (4) 

where 

Pd = PIT -^7(3^-2^ + 3) (5) 

and n s denotes the number of strange quarks in the diquark. 

An approach respecting SU(6) flavour-spin symmetry instead is currently 

investigated. 



5 Our treatment of diquarks resembles to some extent the one employed in the 
Lund approach for baryon production [28]. 
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4 Cluster transitions into primary hadrons 



Once the clusters have been formed, their masses are distributed continuously 
and independently of the hard process with a peak at low mass. In contrast, 
the observable hadrons have a discrete mass spectrum and, hence, the clusters 
must be converted. This is achieved through binary cluster decays and through 
transformations of individual clusters into single primary hadrons. Our model 
does not incorporate the subsequent decays of unstable hadrons. To model the 
cluster transitions, the following assumptions are employed: 

(1) Cluster fragmentation is universal, i.e. independent of the hard process 
and of the parton shower. The clusters disintegrate locally without impact 
on other clusters. 

(2) Cluster transitions, i.e. decays as well as transformations, involve only 
low momentum transfer, of the order of 1 GeV [14], since hadronization 
effects are supposed to be sufficiently soft and event-shape variables such 
as the thrust scale inversely with the centre-of-mass energy. 

(3) The regime of clusters is separated from the regime of hadrons accord- 
ing to the flavours of the cluster constituents and the accessible hadron 
masses. Clusters are supposed to be hadrons, if their mass is below a 
threshold mass. This bound is given by the maximum of the heaviest 
hadron with identical flavour content and the sum of the masses of the 
lightest possible hadron pair emerging in the decay of those clusters. 

The last assumption has two consequences, namely that in a first step the 
newly formed clusters that are already in the hadronic regime have to be 
transformed into hadrons; in the subsequent binary decays of the remaining 
clusters, that also the daughter clusters, which fall into the regime of hadron 
resonances, have to become hadrons immediately. 

In both definite hadron species Ti has to be chosen according to 

the flavour structure of the considered cluster C. Respecting fixed particle 
properties, this choice is based on hadron wave functions motivated by a non- 
relativistic quark model. The wave functions are factorized into a flavour- 
and a spin- dependent part. In our model the flavour part is given for a two- 
component system in terms of quarks and diquarks. The overlap of this flavour 
part with the flavour content of the cluster gives rise to a flavour weight. 
In addition, since spin information is washed out in the clusters [14], the 
total spin J of the hadron manifests itself as a corresponding weight. The 
total spin is given through the coupling of the relative orbital momentum 
L with the net spin S of the valence components. This can be written as 
J = L + S . The contributions of states with different orbital momentum L to 
the total-spin sum are accounted for by some a-priori weights Vl, which enter 
as model parameters. Taken together, the total flavour-spin weight W for a 
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single hadron reads 

W(<?iSi, q 2 s 2 -> 7i J {qiSi, q 2 s 2 ) 





(^F(^ j )|gisi,g 2 s 2 ) 




(<S , |sis 2 ) 


2 


(J\LS) 


2 




(^F(H J )|giSi,g 2 s 2 ) 


2 f 


|sis 2 ) 


2 


(J\LS) 2 



(6) 



In contrast to q\ denoting the quarks, q 2 stands for antiquaries as well as 
diquarks. The spins of the two cluster components 1 and 2 are given by s\ = 
s(qi) and s 2 = s(q 2 ), respectively, and (iPf(H j )\ denotes the flavour part of 
the hadron wave-function 6 . Moreover, |(j|/s)| 2 = 2j + l/(J2i = \i- s \,. ..Ji+s) 
and J2' l s->j * s an abbreviation denoting a summation over L = 0, 1, . . . and 
S = \si — s 2 \, . . . , (si + s 2 ), considering the condition that only those terms 
contribute, where \L — S\ < J < (L + S) can be fulfilled. Finally, it should be 
stressed that the second term of Eq. (6) represents only a static model, which 
accounts for the correct selection of hadrons according to their total spin. 
The cluster fragmentation into primary hadrons is performed in two phases: 
(I) When the clusters are formed from colour-connected pairs of quarks 
and diquarks, some of them, because of their comparably low mass, fall into 
the hadronic regime. Within our framework these clusters are transformed 
into single hadrons immediately. In doing so, however, some four-momentum 
is released and has to be absorbed by other clusters. By allowing hadrons with 
masses lower than the cluster mass only, the momentum transfer is taken 
to be a mere energy transfer and, therefore, is time-like. This ensures that 
the absorbing cluster becomes heavier. To fulfil the low momentum-transfer 
requirement, the already outlined hadron-selection procedure according to the 
flavour-spin weights W is modified through the inclusion of an additional - 
kinematic - weight, which behaves like 



kin. 



exp 



(7) 



In this equation Q 2 > denotes the squared momentum (i.e. energy) transfer, 
and Qo is the scale related to the low momentum-transfer demand. The limit 
Qo, furthermore, depends on the cluster mass and is also employed in the clus- 
ter decays; see below. Note that in the Webber model the clusters being too 
light to decay are identified to be the lightest hadron with identical flavour 
structure [24]. In comparison with the Webber scheme, the major difference 
of our approach in the case of single-cluster transitions is the expansion of the 
hadron-selection procedure. 

The cluster compensating the residual four-momentum is selected such that 
it contains the partner that emerged in the same non-perturbative gluon- 
splitting process as one of the constituents of the transformed cluster. In turn, 



6 For mesons this also includes the possibility of singlet-octet mixing occurring in 
hadron multiplets. 
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Fig. 2. Direct and crossed flavour arrangement and colour flow guaranteeing colour 
neutrality for each final-state configuration in cluster two-body decays. 

clusters, which fall into the hadron regime and contain two leading quarks, are 
always split non-perturbatively into two clusters containing only one leading 
constituent. In this context leading partons, however, are only those quarks 
and antiquarks that directly originate from the perturbative phase, and not 
from the non-perturbative gluon splitting or from the cluster decays. For the 
resulting single-leading clusters, then, the same considerations as for the direct 
transformation to hadrons apply. In case a cluster in the hadron regime is made 
of a diquark and an antidiquark, which is, in principle, possible, it is forced to 
specifically decay into two mesons. The details on the forced double-leading 
cluster breakup and the double-diquark cluster decay are outlined below in 
paragraph (II); see Eq. (8). 

(II) Finally all remaining primary and secondary (daughter) clusters have 
to be split. The mass categorization outlined above automatically yields one 
of the modes C — > CiC 2 , C —> Ci7Y 2 , C — > 7YiC 2 , or C — > Ti.1H.2- These modes 
involve the creation of an extra flavour pair according to the ideas illustrated 
in Sec. 3. Similarly to the cluster-formation phase, then, two flavour configura- 
tions for the decay products emerge, namely a direct one and a crossed one; see 
Fig. 2. Again, the crossed configuration is suppressed by the colour factor 1/Nq 
and the kinematical weight from Eq. (1) using identical measure functions w 
and replacing to by Qq, which again depends on the mass of the decaying 
cluster. The cluster-decay kinematics, which makes use of the parameter Q , 
is fixed to be anisotropic. Starting from a mother cluster with constituent 
momenta 2 , the new momenta of the decay-products' constituents read [14] 



where / and / label the momenta of the newly created flavour pair. Hence, 
for the two cluster arrangements (see Fig. 2) the momenta are given by = 
Pi + Pf, Plr. = Pf + P2 in the direct, and P* oss = p 1 + p 2 , P^ oss = p f + pj in 
the crossed case, respectively. To guarantee well-behaved four-momenta in this 
fission breaking, our model uses a running Qo depending on two parameters, 
Qo and Mo, with the constraint Qo < Mo : 




Qo c 





Qo(M c ) 



Qo-Mc 
M + M c 



< M c . 



(9) 
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Table 1 

Different cluster types emerging through cluster breakups. Decay channels indicat- 
ing a four-quark, i.e. two-diquark, system to become a hadron within the modes 
C — > C1H.2, H1C2, 'H1H2 are vetoed. The four-quark cluster disintegration into two 
mesons (see the last row of the table) is only available for the mode C — > Ti.1Ti.2- 
The occurrence of the two disintegration possibilities is taken to be equally likely. 



Cluster 




Direct case 


Crossed case 




Direct case 


Crossed case 


qi Q2 


qq ^ 


qiq + qg~2, 


qiq 2 + qq 


d75 


qi D + Dq 2 , 


qiq 2 + DD 


qiD 2 


m 


qiq + qD 2 , 


q\D 2 + qq 


DP 


qi D + DD 2 , 


qi D 2 + DD 


D X D 2 


qq 


D 1 q + qD 2 , 


D\D 2 + qq 


DP 


DiD + D~D 2 , 


D\D 2 + DD 


D X D 2 




?2?i + q' 2 q\ 






q2q[ + q' 2 qi 





Having fixed the primary kinematics, via Eq. (8), and the combination of 
flavours and momenta to the new clusters, their masses can be deduced from 
the squares of their total four-momenta. Then, as stated above, the different 
decay modes C — > CiC 2 , C\TC 2 , T~iiC 2 , H.1H2 are distinguished. All possible 
decay channels within each mode are comprehensively summarized in Table 1. 



(1) For the case of breakups involving clusters only, i.e. for C — > C1C2, nothing 
has to be done in addition. 

(2) If one of the daughter clusters falls into the hadronic regime, i.e. for 
C — > C{H.2 and C — > H1C2, a suitable hadron has to be selected such 
that the hadron will be lighter than the cluster. The selection procedure 
follows the one outlined above for the C — ► Ti, transformation; the recoil 
is taken by the daughter system, which belongs to the cluster regime. 

(3) If both new clusters fall into the hadron regime, i.e. for purely hadronic 
decays C — ¥ H1H2, more severe manipulations are applied. First of all, 
the newly created flavour pair // is abandoned; instead, two hadrons 
are chosen directly. Then the combined weight for the selection of such 
a hadron pair consists of three pieces. The first part accounts for the 
two flavour-spin contents. The second one includes the correct relation 
of direct to crossed decay configurations and, furthermore, represents the 
incorporation of the pair-production rates. The last part considers the 
phase space of the decay, which is taken to be isotropic in the cluster's 
rest frame [12,14]. The combination of the first two weights for the hadron 
pair is set up as if only complete SU(3)p multiplets were accessible. Be- 
cause of the superposition with the phase-space factor, a hadron pair that 
cannot be produced in a cluster decay owing to its large mass cannot con- 
tribute to the selection 7 . The other manipulation, as indicated above, is 
that once the hadron species are chosen, the cluster decays isotropically 

The weight treatment for hadron selection in HERWIG was modified by Kupco [29] . 
However, currently, the HERWIG++ group is working on a new approach [18]. 
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in its rest frame into these hadrons. 

Two comments are in order here: first of all, our approach takes leading- 
particle effects into account in the same manner as in Webber's model 
[24]. Secondly, when considering a cluster consisting of two diquarks, 
mesons can emerge only by recombining the individual quarks and anti- 
quarks that constitute the diquarks, see Table 1. Since baryons appear in 
a decay of such clusters through the creation of a quark pair, the diquark 
recombination is taken to be suppressed by a factor of w.r.t. the baryon 
production, which appears with 1—pe in this channel. The specific order- 
ing of the quarks into mesons is then done in a fashion similar to the one 
above, involving hadron pairs. The difference, however, lies in the fact, 
that Clebsch-Gordan coefficients are additionally employed. These coeffi- 
cients account for the rearrangement of the diquark spin wave-functions 
into a double-mesonic basis. 



5 Preliminary results 

The performance of the model introduced above is now illustrated by present- 
ing some results for e + e~ annihilation at the Z° pole using only light quarks 
throughout the event's evolution. The outcomes have been obtained with the 
parton shower of APACIC++-1.0 [30], the matrix elements are generated by 
AMEGIC++-1.0 [31] and matched with this parton shower [32], the primary 
hadronization is accomplished by the cluster model described above, and the 
hadron decays are provided through interfacing the corresponding routines of 
PYTHIA-6.1 [33]. The resulting event generator is the combination of these 
modules. In the following it is referred to as SHERPAa. All results shown below 
are achieved with the same parameter set, where the cluster- model parameters 
have been adjusted by hand. The settings of the other module's input vari- 
ables are mainly taken over from a tuning of APACIC++-1 . 0, together with the 
full hadronization of PYTHIA-6.1. Since measurements that specifically con- 
centrate on the observation of light-quark characteristics are rarely available, 
our results are mainly compared with those gained from running PYTHIA-6 . 1 
and HERWIG-6 . 1 [34] both restricted to u, d, s quarks. Thereby either of the 
models has been run with its default parameter values. 

To begin with, the effects of our colour-reconnection model on the cluster-mass 
distribution, and the statistics of the reconnections in the cluster formation 
are briefly discussed. Figure 3 illustrates the statement that under the influ- 
ence of colour reconnection our cluster hadronization tends to produce less 
massive primary clusters than without the reconnection procedure. The de- 
crease is especially caused by the kinematical factor, Eq. (1), where Wij = p±_^ 
has been used. In the cluster formation one gets approximately 0.74 recon- 
nections per event and, with a frequency of 48%, 35%, 13%, 3%, and 1%, one 
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Primary cluster mass distribution 




M c 



Fig. 3. Primary cluster-mass distribution in electron-positron annihilation events 
that evolve into light-quark and gluon jets at the Z° pole. The SHERPAa result is 
shown with (solid line) and without (dashed line) colour-reconnection (CR) model. 

finds 0, 1,2,3, and > 3 exchange(s), respectively. Moreover, switching off the 
colour-reconnection option entirely while keeping all the other parameters un- 
changed yields the following qualitative modifications: the number of daughter 
clusters per event is increased, which results in an enlargement of the mean 
charged-particle multiplicity of roughly 0.2 tracks per event. The charged-pion 
production rate increases whereas the charged-kaon rate and the (anti)proton 
rate decrease. Furthermore, the charged particle transverse-momentum distri- 
butions are lowered for high p l ^ ont ; however, for the scaled- momentum distri- 
bution of charged tracks, see also below, the bump at x^ ds ~ 0.5 enhances and 
its tail tends to become harder. 

The overall charged-particle multiplicity distribution is presented in Fig. 4. 
The shift to higher multiplicities of the SHERPAa curve w.r.t. the other curves 
indicates the higher mean value of the SHERPAa prediction. Table 2 shows 
mean multiplicities (M^ s ) as provided by those three fragmentation models 
in comparison with inclusive measurements. To exemplify the charged hadron 
rates, the mean multiplicities for the eventually observable charged hadrons 
- n^, and p,p - are considered and compared with experimental uds re- 
sults; see also Table 2. In view of these comparisons, one can conclude that the 
obtained SHERPAa multiplicity results are satisfactory, and in good agreement 
with the PYTHIA-6. l(uds) predictions as well as with the data. 
As an example for a particle-momentum distribution the scaled momentum 
x uds _ 2\p uds \/E cm and its negative logarithm ^ ds = — ha.Xp ds are considered. 
The Xp ds distribution obtained with SHERPAa is shown in Fig. 5, and com- 
pared with the predictions of the PYTHIA-6 . 1 (uds) and HERWIG-6 . 1 (uds) 
event generators. Furthermore, experimental results delivered by the OPAL, 
DELPHI and SLD collaborations on this differential cross section are included. 
The PYTHIA-6 . 1 (uds) model is the most consistent with the OPAL and DEL- 
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Charged particle multiplicity distribution 
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Fig. 4. Predicted multiplicity distribution of charged particles in e + e~ annihilation 
for light-quark and gluon jets at the Z° pole. The SHERPAa result is compared with 
the default PYTHIA-6 . 1 (uds) and HERWIG-6 . 1 (uds) predictions. 

PHI data, but it predicts a slightly softer spectrum. Both cluster-hadronization 
models show a similar behaviour concerning their deviation from these data. 
For Xp ds < 0.7 they wiggle around the PYTHIA-6 . 1 (uds) prediction and for 
x uds > Q g t ne y anticipate a steeper decline, which is quite different from that 
seen in the OPAL and DELPHI data. However, recently published SLD re- 
sults on this topic support the tendency of having a much weaker high-a^ 5 

Table 2 

Overall mean charged-particle multiplicity, and production rates of charged pions, 
charged kaons and (anti)protons in e + e~ collisions. The values are taken for uds 
events running at the Z°-peak centre-of-mass energy. The errors indicated in the 
table are the total errors of the measurements. The abbreviations PY61(nds) and 
HW61 (uds) stand for PYTHIA-6 . 1 (uds) and HERWIG-6 . 1 (uds) , respectively. More 
JETSET and HERWIG results on the topic can be found in [35]. 









(A® 




PY61 (utis) 


19.84 


16.72 


2.010 


0.856 


HW61 (uds) 


18.86 


15.37 


1.693 


1.568 


SHERPAa 


20.15 


16.83 


2.018 


1.047 


OPAL [36] 


20.25 ± 0.39 








DELPHI [37] 


20.35 ±0.19 








DELPHI [35] 


19.94 ±0.34 


16.84 ± 0.87 


2.02 ±0.07 


1.07 ±0.05 


SLD [38] 


20.21 ±0.24 








SLD [39] 


20.048 ± 0.316 


16.579 ± 0.304 


2.000 ± 0.068 


1.094 ± 0.043 
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Charged particle scaled momentum distribution 




Fig. 5. Scaled momentum distribution of charged particles for E cm = 91.2 GeV in 
e + e~ annihilation considering only the light-quark sector. The SHERPAa prediction is 
compared with experimental light-quark data provided by the OPAL, DELPHI and 
SLD collaborations, and to the PYTHIA-6 . 1 (uds) and HERWIG-6 . 1 (uds) outcomes, 
using their default settings. Concerning the mean value (Xp ds ) of the distributions, 
only the HERWIG-6 . 1 (uds) prediction is consistent with the OPAL measurement of 
(x% ds ) = 0.0630 ± 0.0003(stat.) ± 0.0011(syst.) given in [36]. 

tail. This behaviour then is well described by our cluster model. The hump 
at Xp ds ~ 0.5 is truly a deficiency of cluster approaches. In comparison with 
the HERWIG-6 . 1 (uds) prediction, our model yields a smaller bump, and the 
values for Xp ds > 0.9 do not fall off as rapidly as the HERWIG-6. l(uds) ones. 
This performance might be due to the mass categorization treatment of the 
cluster transitions, which has been introduced in our model. All in all, the Xp ds 
behaviour clearly reflects two cluster-model weaknesses, namely (1) that the 
necessary increase in cluster and, therefore, in hadron multiplicity excessively 
results in a decrease of large three-momenta of primary clusters, and (2) that 
the hadronization of events with a small number of primary clusters is not 
sufficiently modelled yet. 

In contrast to the Xp ds distribution, the distribution emphasizes the soft 
momenta of the spectrum. Figure 6 illuminates the SHERPAa result together 
with those of the other two QCD Monte Carlo models, and compares them 
with experimental measurements from the OPAL, DELPHI and SLD collabo- 
rations. SHERPAa describes the data over most of the C,p ds region, and is in quite 
good agreement with the PYTHIA-6 . 1 (uds) prediction. It underestimates the 
region of 1 < £p ds < 2; on the other hand, however, it slightly overestimates the 
data for 3 < £p ds < 5. Experimental inclusive measurements of the peak posi- 
tion, Q uds = 3.76 ±0.02 (DELPHI [35]) and Q uds = 3.74 ±0.22 (OPAL [36]), 
seem to be reproduced by the PYTHIA-6 . 1 (uds) and SHERPAa Monte Carlo 
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Charged particle log( 1 / x ) distribution 




e + e _ annihilation, considering the light-quark sector only. The SHERPAa prediction 
is presented together with experimental uds data provided by the OPAL collabora- 
tion, and with results from default PYTHIA-6 . 1 (uds) and default HERWIG-6 . 1 (uds) . 



simulations. HERWIG-6 . 1 (uds) is considerably low (high) for 2 < £~ s < 5 
(0.5 < ^ ds < 1). 

As an example for the group of event-shape observables, the 1 — T distribu- 
tion, T being the thrust, of the three aforementioned QCD Monte Carlo event 
generators with u, d, s quark restriction is presented in Fig. 7 for light-quark 
and gluon jets. HERWIG-6 . 1 (uds) accounts on average for more spherical event 



Charged particle 1 -thrust distribution 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 
^— mean: 0.0664442 mean: 0.0658546 mean: 0.0695405 -|-T uds 



Fig. 7. 1 — T/xp distribution of charged particles for E cm = 91.2 GeV in e + e~ 
annihilation with a restriction on u, d, s and gluon jets. The SHERPAa prediction is 
compared with predictions of default PYTHIA-6 . 1 (uds) and HERWIG-6 . 1 (uds) . 



14 



Charged particle Durham differential 3-to-2 jet rate 



T3 



10 
10 
1 

10* 
10 
10 
10 




10 



cm energy = 91.2 GeV 



uds only 



MC simulations 

^— PYTHIA-6 . 1 
HERWIG-6 . 1 
SHERPA preliminary 



_L 




0.05 0.1 0.15 0.2 

mean : 0.01 9683 mean : 0.01 92774 



0.25 0.3 0.35 0.' 
mean: 0.0204074 y cu t 



Fig. 8. The Durham 3 — ► 2 differential jet rate of charged particles in e + e~ annihila- 
tion at the Z° pole. Only uds events are taken into account. The SHERPAa result is 
compared with the results stemming from PYTHIA-6 . 1 (uds) and HERWIG-6 . 1 (uds) 
performances, both of which run with their default parameters. 

shapes, which is indicated by a weaker decline of the spectrum towards higher 
values. Owing to the LPHD concept, the SHERPAa prediction, somewhat ex- 
ceeding the PYTHIA-6 . 1 (uds) result for 0.1 < 1—T < 0.3, rather resembles the 
prediction of PYTHIA-6 . 1 (uds) , which might be due to the fact that SHERPAct 
employs a PYTHIA-like parton shower. 

Lastly the Durham 3^2 differential jet rate is considered in Fig. 8. Except 
for the low-statistics region, the results for the event generators shown in the 
plot barely exhibit any deviation from one another. 

Taken together, one can conclude that a good performance could be achieved 
by the new model in electron-positron annihilation into light-quark and gluon 
jets at the Z° pole, although this model has not been sufficiently tuned yet. 
These first SHERPAa outcomes indicate an encouraging agreement with re- 
sults obtained from PYTHIA-6 . 1 restricted onto the light-quark sector. Where 
provided, the comparison with experimental data is satisfactory. 



6 Summary and conclusions 



A modified cluster-hadronization model has been presented. In comparison 
with the long-standing Webber model, the extensions of our approach are the 
following. 

Soft colour-reconnection effects are included in the cluster formation as well as 
in the cluster-decay processes. This yields an enhancement of the number of 
decay configurations. The spin of diquarks is explicitly accounted for through- 
out the model. The number of basic cluster species is enlarged, especially by 
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a new mesonic-cluster type, the four-quark cluster. The significant feature of 
our approach is the flavour- dependent separation of the cluster and hadron 
regimes in terms of the mother cluster's mass. This categorization automati- 
cally selects the cluster-transition mode. Taken together, these aspects require 
the set-up of generically new cluster decay channels. 

Our cluster-hadronization model is implemented as a C++ code. The result- 
ing version is capable of describing electron-positron annihilation e + e~ — > 
7*/Z° — > dd, uu, ss into light-quark and gluon jets. Some first tests were passed 
(see previous section) and the agreement with PYTHIA-6 . 1 (uds) and exper- 
imental data is satisfactory. Some cluster-model shortcomings, such as the 
too low charged-particle multiplicity, could be cured; and the spectrum of the 
scaled momentum could be improved. The model will soon be completed by 
including heavy-quark hadronization. Furthermore, the focus of future work 
is on treating the fragmentation of remnants of incoming hadrons, especially 
in view of proton-(anti)proton applications (Tevatron and LHC physics). 
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